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Ground-based lidars are suitable for long-term ozone monitoring as a complement to satellite and ozo- 
nesonde measurements. However, current ground-based lidars are unable to consistently measure ozone 
below 500 m above ground level (AGL) due to both engineering issues and high retrieval sensitivity to 
various measurement errors. In this paper, we present our instrument design, retrieval techniques, and 
preliminary results that focus on the high-temporal profiling of ozone within the atmospheric boundary 
layer (ABL) achieved by the addition of an inexpensive and compact mini-receiver to the previous system. 
For the first time, to the best of our knowledge, the lowest, consistently achievable observation height has 
been extended down to 125 m AGL for a ground-based ozone lidar system. Both the analysis and pre- 
liminary measurements demonstrate that this lidar measures ozone with a precision generally better 
than ±10% at a temporal resolution of 10 min and a vertical resolution from 150 m at the bottom of 
the ABL to 550 m at the top. A measurement example from summertime shows that inhomogeneous 
ozone aloft was affected by both surface emissions and the evolution of ABL structures. © 2013 Optical 
Society of America 

OCIS codes: (280.1120) Air pollution monitoring; (010.1280) Atmospheric composition; (280.1910) 

DIAL, differential absorption lidar; (010.3640) Lidar; (010.4950) Ozone. 

http://dx.doi.org/10.1364/AO.52.003557 


1. Introduction 

Vertical profiling of atmospheric boundary layer 
(ABL) ozone is important for understanding the inter- 
action between surface and upper-altitude air, 
validating satellite retrievals, and evaluating photo- 
chemical or air-quality models. Anthropogenic 
pollution generally originates at ground level and 
subsequently impacts the downwind surface air 
quality through transport over various spatiotempo- 
ral scales [1-3]. Surface ozone is coupled with the 
upper-air ozone, although these regions have distinct 
diurnal variations. 
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As a complement to satellite and ozonesonde mea- 
surements, lidar profiling in the ABL has been pre- 
viously carried out using both airborne [4-6] and 
ground-based systems [7—9]. These lidars generally 
retrieve ozone profiles with a vertical resolution of 
hundreds of meters and an integration time of mi- 
nutes. These efforts have significantly added to our 
understanding of tropospheric chemistry and dy- 
namics. However, the quantity of these ozone lidar 
data has been limited and unable to satisfy our re- 
quirement for high-resolution data, especially within 
the ABL where both the vertical resolution of the 
satellite retrieval and the temporal resolution of 
the ozonesonde are coarse. 

Although airborne lidars retrieve near-surface 
ozone over large spatial regions, they have limited 
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precision near the surface, which is at their far range. 
Moreover, airborne lidars require significant opera- 
tional costs that limit their deployment. To our 
knowledge, except for lidars operating in a horizontal 
scanning configuration [10], the lower limit for ozone 
retrievals by ground-based lidars is generally 500 m 
above ground level (AGL) due to several difficulties. 
First, near-surface ozone measurement requires 
rapid overlap of the receiver field of view (FOV) 
and laser beam divergence. This overlap can signifi- 
cantly increase the near-range signal intensity that 
saturates the detector and results in nonlinear sig- 
nals [ 11 ] . The conventional chopper or gating system 
response is not rapid enough to protect the system. 
Second, the large dynamic range, due to strong 
Rayleigh and aerosol scattering, requires multiple 
detectors to cover the entire ABL (< ~ 2 km). This 
adds to the complexity of the detector design and 
the retrieval algorithm. Third, synchronization of 
the timing between on- and off-line wavelengths is 
critical because of the rapidly decaying backscatter 
signal. Small timing offsets will result in large ozone 
retrieval errors. Fourth, the uncertainty in the 
lower ABL ozone lidar retrieval caused by aerosol 
interference can be high due to differential aerosol 
backscattering, and the unknown Angstrom expo- 
nent and lidar ratio. 

In this article, we point out the challenges for 
ozone lidar measurement within the ABL, particu- 
larly near the surface, and present our design 
solutions and retrieval methodology for overcoming 


these difficulties. Then, we show an example of our 
ozone measurements by focusing on the ABL. Fi- 
nally, we examine the resulting retrieval precision. 

2. Instrument 

Our tropospheric ozone differential absorption lidar 
(DIAL) is located on the fourth floor of the National 
Space Science & Technology Center (NSSTC) on the 
campus of the University of Alabama in Huntsville 
(UAHuntsville) at an elevation of ~206 m above 
sea level. This system currently retrieves ozone be- 
tween 0.125 and ~12 km AGL using three receivers 
(2.5, 10, and 40 cm) and a wavelength pair of 285 
and 291 nm. 

A. Laser Transmitter 

Table 1 presents the characteristics of the UAHunts- 
ville ozone lidar. The laser transmitter consists of two 
wavelength-tunable dye lasers, pumped by two sepa- 
rate frequency-doubled Nd:YAG lasers at 532 nm. The 
30 Hz pump lasers are triggered by an external pulse 
generator with a 16.7 ms separation between the al- 
ternate pulses. The dye lasers are frequency doubled 
to generate two ultraviolet (UV) laser beams at 285 
and 291 nm with energy of 5-7 mJ • pulse -1 and pulse 
length of ~7 ns. The outgoing beams have near-field 
diameters of ~6 mm and divergences of less than 
1 mrad. Compared with the previous system [12], 
the pump lasers were upgraded in 2010 and the 
resulting increase in UV energies raised the upper 
observation limit to ~12 km. 


Table 1. Characteristics of the UAHuntsville Ozone Lidar 


Transmitter 


Specification 


Pump lasers 

Nd:YAG, 30 Hz repetition rate, 150 mJ • pulse 1 at 532 nm 


Dye 


Mixture of rhodamine 590 and 610 


Emitted UV lasers 

7 mJ ■ pulse 1 at 285 nm, 5-6 mJ • pulse 1 

at 291 nm; 0.6 cm near-field beam diameter, 


7 ns 

pulse length and divergence <1 mrad for both lasers 


Receiver 

Channel 1 

Channel 2 

Channel 3 

Channel 4 

Diameter (cm) 

2.5 

10 

10 

40 

Focal length (m) 

0.1 

2.3 

2.3 

4.5 

Separation from the laser beams (cm) 

10 

20 

20 

50 

FOV (mrad) 

10 

4.3 

4.3 

1.5 

Full overlap height (meter above 

20 

90 

90 

1200 

the laboratory) 





Light split percentage (%) 

100 

10 

90 

100 

Solar blind filter 

Center wavelength at 286.4 

nm with all nm FWHM; 35% 



transmittance at 285 nm and 20% transmittance at 291 nm 

Neutral density filter transmittance 

10% 

10% 

32% 

NA 

PMT type 

Hamamatsu R9880U 

Hamamatsu R7400U 

Hamamatsu R7400U Electron Tube9831QA 

Gated PMT delay (ps) 

0 

0 

1 

>10 

Measurable height range (km) 

0.1-1 

0.4-1. 5 

1-4 

3-12 

Signal processing 

Channel 1 

Channel 2 

Channel 3 

Channel 4 

Photon counting 


250 MHz maximum counting rate 


Analog 


12 bit and 40 MHz analog-to-digital converter 


Fundamental range resolution (m) 

3.75 

3.75 

3.75 

3.75 

Degraded range resolution 

22.5 

45 

75 

150 

due to average (m)“ 






“This resolution is equal to the range interval to report ozone. In terms of the effective spatial resolution of the DIAL retrieval, please 
refer to Fig. 5. 
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The selection of on- and off-line wavelengths at 285 
and 291 nm is a compromise between minimizing the 
solar background, reducing the impact of the aerosol 
interference, and maximizing the range limit. 
Although the uncertainties in the DIAL retrieval 
due to aerosols have a complicated relationship with 
the wavelengths, their separation, and the local aero- 
sol conditions, wavelengths between 285 and 300 nm 
have been previously shown [13-16] to retrieve ozone 
with acceptable errors if the aerosol correction is 
treated appropriately [5,17-19]. 

B. Lidar Receiver 

The receiving system consists of three receivers and 
four photomultipliers (PMTs) as shown in Fig. 1 and 
listed in Table 1. The 2.5, 10, and 40 cm receivers 
fully overlap both laser beams at 20, 90, and 
1200 m, respectively, above the laboratory height. 
The lidar receiver configuration is similar to that re- 
ported by Kuang et al. [12] except for the addition of a 
2.5 cm mini-receiver and a beamsplitter in the 10 cm 
receiver. In 2008, a 10%-90% beamsplitter was inte- 
grated into the 10 cm receiver and extended the ob- 
servations down to ~500 m AGL using the 10% 
channel [20,21]. Modifying a design used in a pre- 
vious National Aeronautics and Space Administra- 
tion (NASA) airborne lidar campaign [22], a 2.5 cm 
diameter receiver was recently added to lower the li- 
dar’s minimum observable altitude. This approach, 
using four PMT (or altitude) channels, provides nec- 
essary dynamic range to cover both the ABL and 
the free troposphere under expected atmospheric 




Fig. 1. Configuration schematic for the UAHuntsville ozone lidar 
in (a) its vertical cross section and (b) its horizontal cross section. 


conditions. For each laser, a photodiode detects the 
outgoing laser pulse and triggers a function genera- 
tor to control the timing sequence for both the Licel 
(Licel GmbH, Germany) transient recorder (TR, 
model 40-80) and the gating system. The PMTs 
for the two high-altitude channels (channels 3 and 
4) are gated typically for the first 1 and 10 ps, respec- 
tively, to avoid saturation. Although this paper 
discusses only the ABL ozone measurement, we still 
present the parameters for all of the channels in 
Table 1 for convenient comparison. 

As shown in Fig. 2, the mini -receiver consists of an 
adjustable focusing element holding a 10 cm focal 
length plano-convex lens that facilitates the optimi- 
zation of the distance between the lens and the PMT. 
An adjustable iris is inserted between the 10 cm lens 
tube and PMT to control the receiver’s FOV. The iris 
can be replaced by a suitably sized pin hole after the 
FOV is experimentally determined. Currently, the 
FOV is set to 10 mrad to ensure that the receiver 
and laser are fully overlapped at ~20 m. At the lower 
part of the mini-receiver, a Licel module incorporat- 
ing a compact Hamamatsu R9880U PMT detects 
backscattered photons. Although the photocathode 
area size of this R9880U PMT is only 8 mm, it is suf- 
ficient to fully cover the backscattered light that 
passes through the receiver aperture without a diffi- 
cult manual alignment. A solar filter is employed for 
daytime measurements. A 10 dB neutral density fil- 
ter is typically added into channel 1 to maintain the 
linearity of backscatter signals because its PMT is 
not gated. 

C. Range Registration Offset 

For ozone retrievals, the ranges of the on- and off-line 
wavelengths should have the same reference height, 
which is typically the ground. The offset in range 
registration between the two wavelengths, r* (or t* = 
2 r*/c accounting for the two-way travel, where c 
is the light speed), was not emphasized for high- 
altitude retrievals in previous publications because 
of the slowly decaying signal and long integration 
range interval (usually more than 100 m). However, 
minimizing this offset is critical for near-surface 
measurements. 



Fig. 2. Design of the 2.5 cm receiver. 
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Fig. 3. Timing chart for the 2.5 cm receiver. r PD represents 
the range (or the height AGL) of the photodiode trigger pulse, 
and r s denotes range at the start of the data recording by the 
Licel TR, which is a photon counter. Subscripts 285 and 291 
represent the two lasers. Receiver is located at approximately 
12 m AGL for both lasers, which is the laboratory height. The solid 
portions of the signals represent the signals that are recorded by 
the TR. The dashed portion is not recorded. For accurate DIAL 
retrieval, r* should be minimized by adjusting Delay 2 s 5 and 
Delay 29 i. 


Figure 3 presents the schematic timing chart for 
our 2.5 cm receiver channel (channel 1). The delay 
(denoted by Delay 285 or Delay 2 gi for the two lasers) 
between the photodiode trigger (also refer to Fig. 1) 
and the start of data recording consists of the time 
that signals take to traverse the cables [typically 
the Bayonet Neill-Concelman (BNC) cables] and 
the internal delay between the different trigger 
channels of the function generator and, finally, the 
internal delay in the Licel TR. Because of these de- 
lays, the TR misses the dashed portion of the back- 
scattered signals in Fig. 3. However, further effort to 
reduce those delays is not necessary because the 
missed signals are below the full overlap altitude 
where we are currently unable to obtain reliable 
retrievals. To minimize r* so that the 285 and 
291 nm signals have the same range registration, 
we adjust Delay 285 and Delay 291 by changing the 
BNC cable lengths and monitoring both Delay 285 
and Delay 291 using an oscilloscope. The pulse lengths 
and shapes observed by the photodiode need to be ac- 
counted for if they are different for the two UV 
beams. In practice, if the pulse shapes are difficult 
to measure, an empirical correction (e.g., tuning 
cable length based on the calibration using other 
instruments) may be necessary for the near-range 
retrievals even after various efforts have been taken 
for minimizing r*. 

Figure 4 shows modeling results for errors in the 
DIAL retrievals due to r* as a function of range by 
assuming a 1.5 x 10 12 molecule cm -3 constant ozone 
profile, no aerosol, and using the 1976 U.S. Standard 
Atmosphere [23]. The ozone number densities were 
computed without any filtering or smoothing because 
no modeled random errors were added onto the simu- 
lated signals. When r* is greater than 0.5 m, the 
retrieval errors become significant (>10%) below 


Errors (%) 



Range offset (m) 

Fig. 4. Modeled percentage errors in the DIAL retrievals due to 
range registration offsets, r * , between on- and off-line wavelength 
by assuming a 1.5 x 10 12 molecule cm -3 constant ozone profile, no 
aerosol, and using the 1976 U.S. Standard Atmosphere [23], Every 
1.5 m range offset corresponds to a 10 ns time offset. 

200 m. The retrieval errors are sensitive to r*, 
especially when r* is small, and insensitive to the 
designated retrieval interval (Ar in DIAL equation). 
These errors are primarily attributable to the differ- 
ent Rayleigh extinctions at ranges between r and 
r + r* and decrease quasi-exponentially with range 
because the Rayleigh extinction decreases exponen- 
tially with range. 

3. Data Processing 

To improve measurement precision, several tech- 
niques are employed on the raw lidar signal. The 
temporal integration time is typically set at 
10 min. The integrated spatial range intervals are 
22.5, 45, and 75 m for PMT channels 1, 2, and 3 
(Fig. 1). The nonparalyzable dead time [ 24 ] has been 
empirically determined to be ~4 ns and is largely as- 
sociated with the type of PMT and the counter’s dis- 
criminator level. Far-range signals, where no laser 
signal returns are expected, are averaged to give 
an appropriate background for each wavelength 
and PMT channel. Although analog signals (12 bit) 
are conventionally considered to be less reliable than 
photon counting (PC), they can increase the virtual 
counting rates and thus improve the dynamic range. 
The analog and PC signals are merged together at 
counting rates between 20 and 30 MHz after their 
offset [25], ~250 ns for our system, is removed. An 
exponential function, determined using a least 
squares fit to data acquired approximately from 
100 to 160 ps, simulates the signal-induced bias 
(SIB) [26], The SIB correction is empirical; however, 
it can significantly impact the precision of the far- 
range retrievals. Two dominate factors affecting 
SIB are the PMT type and the backscattered signal 
intensity. In our measurements, SIB is often easy to 
characterize in the nighttime data when the solar 
background is nonexistent. 
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4. DIAL Retrieval 

The DIAL equation can be simplified as [ 27 ] 


^o3(r) 


-1 

2A0 3w 



P 

1 on(r) 
-^off (r) _ 


+ An o3(r) , 


(1) 


where n,, : > lfr) is the ozone number density at range r, 
A 0 3(r ) is the differential ozone absorption cross sec- 
tion, P 0 n(r) and P 0 ff(D are the backscattered on- and 
off-line lidar returns, and A?i o3() .) is the residual ozone 
term including Rayleigh and aerosol corrections. The 
first term of the right side of Eq. (1) is often called 
the signal term. 

After processing, the ratios of on- and off-line 
signals as a function of range are computed. The 
derivative of the logarithm of the signal ratios can 
be calculated by either the discrete form or an inner 
derivative function of certain advanced computer 
language (e.g., IDL or MATLAB). We choose the 
latter method because it is easier to deal with the 
edge effect of the near-ground signal and is more flex- 
ible to filter out the random noises with variable 
widths of the convolution window. Savitzky-Golay 
(SG) low-pass filters [ 28 ] can smooth the data using 
a polynomial of arbitrary degree and width, and un- 
equal spacing between data points. Therefore, this 
method can handle well the lowest data point and 
is particularly suitable for near-ground signal 
processing. We employ a first-order SG differentia- 
tion filter with a second-degree polynomial and var- 
iable fitting widths to In P on(r) /P 0 \u r) smooth. A large 
width is preferred at far ranges where the signal-to- 
noise ratios (SNRs) are relatively low. Filtering 
widths over the same altitude range can also vary 
with PMT channels, where their detected signal 
intensities are significantly different. 

The effective spatial resolution of the final ozone 
retrieval (A r e ), which determines the ability to see 
fine detail in the data, is not equal to, but is associ- 
ated with the convolution (fitting) window width 
(L s ). Ar e is always better than L s and worse than 
the integrated signal interval (Ar). A r e can be defined 
by Ls/fc [29], where fc is the cutoff frequency at a 
predefined attenuation threshold, typically -3 dB. 
The optimum 3 dB down fc of the SG filter for the 
second-degree polynomial fitting is close to 2 Hz 
[30], theoretically. However, previously reported res- 
olutions in lidar retrievals calculated by various 
methods [31,32] were lower than the optimum value. 
We conservatively assume an fc at 1.6 Hz [33] to es- 
timate the effective resolutions of the DIAL retriev- 
als as a function of range for each channel shown by 
Fig. 5. Generally, a symmetrical (bins on the left side 
and right side are equal) SG filter with a half-width 
of 6 bins is applied on the logarithm of the signal ra- 
tios for all channels. In spite of the same number of 
bins, the filtering widths of the channels still differ 
because channels 1, 2, and 3 have bin lengths of 
22.5, 45, and 75 m, respectively. At the lower-range 
limit where the data points do not allow for a sym- 
metrical window, a filter with a nonsymmetric width 



Fig. 5. Estimated effective spatial resolutions of the DIAL 
retrievals as a function of range for channels 1, 2, and 3. The 
step-function-like increase in resolution at the upper range of each 
channel is due to an increasing filtering width. The gray shaded 
areas designate the region where the adjacent channels can be 
joined together. 


is employed. Because of the higher uncertainties in- 
herent in using a nonsymmetric filter, especially in 
the presence of a large ozone gradient condition, 
the ozone retrievals in the initial 2 or 3 bins are dis- 
carded, and the first valid retrieval occurs at ~122 m 
AGL (or ~110 m above laboratory). At the upper 
range limit a filter with a linearly increasing width 
is applied due to the low SNRs. In the ABL, the esti- 
mated retrieval A r e increases with height from ~150 
to ~550 m. The ozone retrievals of adjacent channels 
can be joined together based on their statistical un- 
certainties [ 12 ] in the designated regions as shown 
in Fig. 5. 

The pressure and temperature profiles for atmos- 
pheric density, Rayleigh scattering, and ozone ab- 
sorption cross section are provided by local 
ozonesonde or radiosonde measurement. Without a 
third wavelength, an iterative aerosol correction pro- 
cedure [ 12 ] is applied to reduce errors arising from 
differential aerosol backscattering. This procedure 
substitutes the first estimated ozone from the signal 
term in Eq. (1) to derive the aerosol profile for the off- 
line wavelength and iterates the aerosol and DIAL 
calculation until obtaining a stable solution. How- 
ever, this procedure needs to assume a value for 
the Angstrom exponent and lidar ratio. These as- 
sumptions result in additional retrieval uncertainty. 

5. Ozone Measurements 

Figure 6 shows a lidar retrieval made on 30th June 
2012 with a 10 min temporal integration (18,000 
shots) and varying spatial resolution. The collocated 
ozonesonde retrieval at 13 : 10 local time and the 
Environmental Protection Agency (EPA) (~16 km 
away from lidar station) hourly surface measure- 
ments are shown as a comparison and a preliminary 
validation. The ozonesondes used in Huntsville 
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Fig. 6. Ozone lidar retrievals (channel 1 — bottom, channel 2 — 
middle, and channel 3 — top) compared with the ozonesonde 
(marked by the black triangle at 13:10 launch time) and EPA 
(~16 km away) hourly surface measurements. 


employ an electrochemical concentration cell sensor 
[34], with 2% unbuffered potassium iodide (KI) 
solution, to measure ozone with a vertical resolution 
of 100,150 m [ 35 ] and an accuracy better than ±10% 
[36,37]. The ozone retrievals from three lidar chan- 
nels (two independent receivers) are generally con- 
sistent with each other over their retrievable 
altitude ranges. All three channels’ retrievals display 
the same ozone reduction below 1 km at 21:00 local 
time. The random noise in the DIAL retrievals in- 
creases with height due to the increasing statistical 
errors, which will be described in Section 6. The 
channel 1 daytime retrievals above 500 m appear 
noisier than nighttime and somewhat inconsistent 
with either channel 2 or the ozonesonde as a result 
of daytime solar background. The solar filter for 
channel 1 has somewhat lower out-of-band blocking. 
Channel 1 also has higher solar-to-signal ratios than 
the other two channels due to a larger FOV. However, 
this does not have a significant impact on the final 
joined retrievals because the altitude above 500 m 
is covered by channel 2. 

A large ozone gradient associated with a temper- 
ature inversion (not shown) on the ozonesonde pro- 
file suggests a mixed layer (ML) with top at 
~2.1 km. The top of the ABL and its evolution can 
also be clearly seen on the lidar curtain plot (Fig. 6), 
indicating that the lidar has captured this important 
feature. Before 18:00, the EPA surface ozone varied 
between 70 and 80 ppbv with a peak at 12 : 00. After 
18 : 00, the surface ozone decreased almost linearly to 
~20 ppbv at midnight because of the development of 
a nocturnal (stable) boundary layer and nitrite titra- 
tion of the ozone. The ozone evolution at higher alti- 
tudes differed from that observed at the surface. The 
upper-air ozone in the ML increased with time and 
peaked at ~90 ppbv several hours later between 
17:00 and 18:00 due to mixing from the surface 



Fig. 7. Comparison of the 10 min DIAL retrievals for the three 
altitude/PMT channels and the collocated ozonesonde measure- 
ment. The mixed layer height is about 2.1 km, suggested by a large 
ozone gradient and confirmed by a temperature inversion (not 
shown) from the ozonesonde profile. 


and in situ production. Surface emissions are the 
principle pollution source typically observed in the 
summer daytime with little horizontal transport. 
However, it is not fully understood why the ozone 
aloft, originating at the surface, has a higher concen- 
tration than at the surface, although this phenome- 
non has been observed previously [38,39] . The ML 
ozone decreased ~ 10 ppbv between 20:00 and 
21:00 mostly as a result of the decreasing photolysis 
rate caused by the intermittent cloud cover above the 
ABL [40]. After 21:00, the ozone below 1 km fell to 
~60 ppbv, while the upper ABL ozone remained at 
80-90 ppbv, suggesting weak mixing in the residual 
layer. 

Figure 7 provides a detailed comparison between 
the 10 min DIAL retrievals for all three channels 
at 13 : 05 and the ozonesonde measurement at 
13 : 10. The ozone number densities on channels 1 
and 2 are consistent in their overlapped region; how- 
ever, they are ~10% lower than the sonde at some 
heights. The channel 1 ozone retrievals agree with 
the sonde below 500 m and are slightly lower than 
the sonde above 500 m. Slight adjustment of the 
SIB correction can correct the channel 1 retrievals 
at the upper range for this particular time. However, 
based on the whole lidar curtain in Fig. 6, those 
biases appear to be due to random noise originating 
with the uncharacterized solar background as men- 
tioned previously. Interestingly, both the ozonesonde 
and lidar show a large ozone gradient (more obvious 
when displayed in units of mixing ratio) below 
~300 m in spite of a sunny summer midday. These 
measurements demonstrate that the upper-air ozone 
variations can be significantly different from the sur- 
face and the ML may not be well mixed. 
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6. Uncertainty Discussions 

Measurement uncertainties can be divided into 
random and systematic errors. The statistical uncer- 
tainty, arising from signal and background noise 
fluctuations, is the principal random error source 
determining measurement precision. Because of its 
importance, the statistical uncertainty has been 
extensively discussed in previous publications 
[ 27,41, 42], Generally, the statistical uncertainties 
increase with height due to the decaying SNR and 
become significant at the far range. With the 
assumption of a Poisson distribution governing PC, 
the statistical uncertainty is approximately propor- 
tional to (ArgAiAP/,) -1 / 2 [ 31 , 42 ], where N represents 
the number of laser shots associated with the inte- 
grated time period, A is the receiver area, and P L is 
the number of emitted laser photons. This means that 
an acceptable retrieval precision can be obtained by 
adjusting the spatial and temporal resolution. 
Figure 8 shows the modeled statistical uncertainties 
for our system with the vertical resolutions in Fig. 5 by 
using the averaged summer ozonesonde profile [ 43 ] . 
The statistical uncertainties are typically less than 
4% within each channel’s measurement range 
(Table 1). 

Systematic errors are associated with (1) differen- 
tial backscattering and extinction due to air mole- 
cules (Rayleigh scattering), aerosols, and other 
gases (e.g., S0 2 , N0 2 ); (2) uncertainties in the ozone 
absorption cross section; and (3) instrumental uncer- 
tainties. The differential Rayleigh scattering results 
in an inaccuracy less than 1% [5] using the pressure 
and temperature profile measured by local ozone- 
sonde or radiosonde. The estimated error due to 
S0 2 or N0 2 for the 285-291 nm pair is minor, less 
than 1% [12]. The aerosol impact can be reduced 
to a minor level for a two-wavelength lidar by cou- 
pling a Klett [44,45] or Fernald [ 46 ] aerosol inversion 
algorithm [5,18,47] into the DIAL retrieval if the 



Statistical uncertainty (%) 


Fig. 8. Modeled statistical uncertainties for channels 1, 2, and 3 
with the effective vertical resolutions shown in Fig. 5 by using an 
average summer ozonesonde profile. 


Angstrom exponent and lidar ratio are exactly 
known. The primary concern arises from the uncer- 
tainties in the assumed Angstrom exponent and 
lidar ratio, which have not been well characterized 
in the UV region (most data is at 308 and 
355 nm), and also the lack of local measurement to 
characterize aerosol types and their properties. 
Kuang et al. [12] have shown that the estimated 
errors due to aerosols in our lidar retrievals are 
within ±(5-20)% depending on the accuracy of the 
Angstrom exponent and lidar ratio. 

Orphal [ 48 ] reviewed several current ozone cross- 
section data sets and found that they have an overall 
agreement to within l%-2% after accounting for 
wavelength shifts, baseline effects, and spectral 
resolution. We choose to use the Brion-Daumont- 
Malicet (BDM) ozone cross-section database [49-52] 
recommended by Liu et al. [ 53 ] for satellite retrievals 
from the UV measurements. The cross-section tem- 
perature dependence, changing slightly (within 
2%) for the Hartley band, is modeled by the quadratic 
polynomial fits from four temperatures 218, 228, 243, 
and 295 K [ 53 ] . Based on Orphal’s [ 48 ] review, we be- 
lieve the uncertainty in the differential ozone cross 
section is less than 2%, translating into a DIAL 
retrieval uncertainty of a similar magnitude accord- 
ing to Eq. (1). 

As shown by Fig. 4, the near-range lidar retrievals 
are sensitive to r*. The retrieval uncertainty due to r* 
can be reduced to less than 1.5% above 125 m AGL 
through use of an oscilloscope with time measure- 
ment precision within 0.2 ns. Moreover, this system- 
atic error does not vary with atmospheric conditions 
or SNR, and thus can be calibrated against local 
ozone observations. We report retrievals only above 
the full overlap altitude for our ozone lidar and do not 
employ any partial overlap function. The SIB impact, 
varying largely with the near-range signal intensity, 
can be significant at the far ranges (>3 km). 
However, it can be reduced to ±5% with empirical 
corrections. 

The total root-mean-square error in the DIAL 
retrieval from all above sources is less than 10% 
excluding an extremely large-gradient aerosol condi- 
tion. For checking the quality of the lidar retrievals, 
Fig. 9 presents the mean and l - a standard 
deviation of the 10 min lidar retrievals in Fig. 6 
for each channel. It is natural that we always prefer 
the retrievals made by the channel with higher pre- 
cision at the overlapped region. The ozone variability 
is associated with both the true atmospheric varia- 
tions and the measurement uncertainties. The stan- 
dard deviations are less than 10% below 4 km except 
for the larger deviations below 1 km mostly due to 
the ozone diurnal variations related to the ABL 
development suggested by Fig. 6. This reasonable 
variability is consistent with the measurement un- 
certainty analysis. 

It is worth noting that the retrieval errors due to 
aerosols are large near inhomogeneous layers, such 
as just above the surface and near the top of the ABL 
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(a) (b) 


Fig. 9. Variability of the ozone lidar measurements in Fig. 6. (a) Average of the 10 min lidar profiles and the 1 — a standard deviation, 
(b) The l-i r standard deviation in units of percentage. 


where the statistical uncertainties are relatively 
small for the ground-based system. Based on the 
above analysis and the comparison of the retrievals 
and the ozonesonde observation, the estimated total 
measurement precision for altitudes up to 4 km AGL 
is about ± 10% under normal aerosol conditions and 
could be up to ±20% under an extremely turbid sky 
condition. This conclusion is consistent with Figs. 7 
and 9, which show that the retrievals from the differ- 
ent channels agree within ±5% and the lidar agrees 
with the ozonesonde within ±10%. 

7. Conclusions 

Previous ground-based ozone lidars that were not 
operating in a scanning configuration were unable 
to consistently measure ozone below ~500 m AGL. 
By adding a 2.5 cm receiver to the previous system 
[12], we have expanded the ozone measurements 
down to 125 m AGL and significantly reduce the 
observational gap between the surface and upper- 
altitude air. We have presented a hardware design, 
signal processing and retrieval techniques, error 
analysis, and a preliminary validation result for 
ground-based ozone lidar measurements within 
the ABL. 

Using two receivers and three PMT (altitude) 
channels, the UAHuntsville lidar measures ABL 
ozone with a typical temporal resolution of 10 min 
and spatial resolutions ranging from ~150 to 
550 m. A careful experiment to minimize the range 
offset (r*) between the on- and off-line wavelength is 
critical for measuring the near-ground ozone. After 
the proper calibration of r*, the lidar was able to 
monitor ozone up to 4 km AGL with errors typically 
within ± 10%, suggested by both the analysis and the 
collocated ozonesonde measurement. 

The lidar measurements in Fig. 6 present a nice 
example of the ABL ozone variations affected by 
the surface emissions and the diurnal variation 
of the ABL structures. The distinct evolutions of 
the surface and upper-air ozone demonstrate the 


importance of the high-resolution ABL ozone profil- 
ing, particularly for air-quality models. 
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